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Abstract 

Efforts to unify the Bardeen, Cooper & Schrieffer (BCS) and the Bose-Einstein condensation (BEC) 
formahsms in terms of a ''complete boson-fermion (BE) model" (CBEM) are surveyed. A vital distinction is 
that Cooper pairs (CPs) are indeed bosons that suffer BEC, in contrast with BCS pairs that are not bosons. 
Another crucial ingredient (particularly in 2D where ordinary BEC does not occur) is the linear dispersion 
relation of "ordinary" CPs, at least in leading order in the center-of-mass momentum (CMM) power-series 
expansion of the CP energy. This arises because CPs propagate not in vacuo but in the Fermi "sea." A 
many-body Bethe-Salpeter equation treatment of CPs based on the ideal Eermi gas (lEG) sea yields the 
familiar negative-energy, two-particle bound-state if 2h-CPs are ignored as in the ordinary CP problem. But 
it gives purely-imaginary energies, and is thus meaningless, if 2h-CPS are included as completness requires. 
However, when based on the BCS ground state instead of the IPG, in addition to the familiar trivial solution 
(or Anderson-Bogoliubov-Higgs) sound mode, legitimate two-particle moving "generalized CPs" emerge but 
as positive-energy, finite-lifetime, resonant nontrivial solutions for nonzero-CMM. This amounts to replacing 
the purely-kinetic-energy unperturbed Hamiltonian by the BCS one. The moving CPs again have a linear 
dispersion leading term. BEC of such pairs may thus occur in exactly 2D (as it cannot with quadratic 
dispersion) and in fact all the way down to (1 -I- e)D where e can be infinitesimally small, thus encompassing 
all empirically known superconductors. 

The unified theory reduces in limiting cases to all the main continuum (as opposed to "spin") statistical 
theories of superconductivity. These include both the BCS and BEC theories. The unified BE theory is 
"complete" in that not only two-electron (2e) but also two-hole (2h) CPs are allowed, and in arbitrary 
proportions. In contrast, BCS theory can be deduced from the CBEM but allows only equal (50%-50%) 
mixtures of them, a fact rarely if ever stressed. The CBEM shows that the BCS condensate is precisely a BE 
condensate of a mixture of kinematically independent electrons coexisting with weakly-coupled zero CMM 
2e- and 2h-CPs in equal proportions. Without abandoning the electron-phonon mechanism, the CBEM has 
been applied in 2D and 3D. The BCS model interaction in moderately weak coupling is sufficient to reproduce 
the unusually high values of Tc (in uruts of the Eermi temperature) of 0.01 — 0.1 empirically exhibited by the 
so-called "exotic" superconductors, including cuprates. This range is high relative to the low values of < 10~^ 
more or less correctly reproduced by BCS theory for conventional (mostly elemental) superconductors. Also 
accounted for is the empirical fact that "hole superconductors" systematically have higher Tc's. Room 
temperature superconductors are predicted to be possible but only via 2h-CP BE condensates. 

Running Title: BCS and BEC Unification: A Review 
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1 Introduction 



Boson-fermion (BF) models of superconductivity (SC) as a Bose-Einstein condensation (BEC) [1] [2] go back 
to the mid-1950's [3]- [6], pre-dating even the BCS-BogoUubov theory [7]- [9]. Although BCS theory only con- 
templates the presence of "Cooper correlations" of single-particle states, BF models [3]-[6],[10]-[18] posit the 
existence of actual bosonic CPs. Such pair charge carriers have been observed in magnetic flux quantization 
experiments on elemental [19. 20] as well as on cuprate [21] superconductors (SCs). Larger clusters than pairs 
are not observed, apparently because the clustering occurs not in vacuo but in the Fermi sea. However, no 
experiment has yet been done, to our knowledge, that distinguishes between electron and hole CPs. CPs ap- 
pear to be the single most important universally accepted ingredient of SC, whether conventional or "exotic" 
and whether of low- or high-transition-temperatures Tc. And yet, inspite of their centrality they are poorly 
understood. The fundamental drawback of early [3]- [6] BF models, which took 2e bosons as analogous to di- 
atomic molecules in a classical atom-molecule gas mixture, is the notorious absence of an electron energy gap 
A(T). "Gapless" models cannot describe the superconducting state at all, although they are useful in locating 
transition temperatures if approached from above, i.e., T > Tc- Even so, we are not aware of any calculations 
with the early BF models attempting to reproduce any empirical Tc values. The gap first began to appear 
in later BF models [10]-[15]. With two [13] [14] exceptions, however, all BF models neglect the effect of hole 
CPs included on an equal footing with electron CPs to give a "complete" BF model (CBFM) consisting of both 
bosonic CP species coexisting with unpaired electrons. 

For 13 years the highest Tc value for any superconductor was 23 K, until the discovery [22] in 1986 of the first 
so-called "high- Tc" cuprate superconductor LaBaCuO having a Tc ~ 35 ii". A feverish search for materials with 
even higher Tc's lead within just seven years to the highest-Tc superconductor known to date, the HgBaCaCuO 
cuprate [23] with a Tc ~ 164 K under very high pressure (~ 310,000 atm). The embarrassing fact is that since 
1993 this record has not been broken, very probably because there is yet no predictive microscopic theory of 
superconductivity that can provide orientation in a search that until now has proceeded on a trial-and-error 
basis. 

We submit that progress in developing such a theory has been held back by too many common myths or 
false dogmas firmly entrenched in the theoretical community. 

2 Some false dogmas 

The false theoretical dogmas just mentioned can be summarized in the following assertions: 

1. With the electron-phonon dynamical mechanism transition temperatures (as given by the BCS formula) 
Tc = 1.136£)e~^/^ < 45 at most, since a typical Debye temperature 0^ 300K and A < ^. For higher Tc's 
one needs to invoke magnons or excitons or plasmons or other electronic mechanisms to provide pairing [24]. 

2. Superconductivity is unrelated to Bose-Einstein condensation (BEC) [25]. 

3. BEC is impossible in 2D [26] [27]. 

4. Cooper pairs: 

a) are such that "there is a very strong preference for singlet, zero- momentum pairs, so strong that one 

can get an adequate description of SC by treating these correlations alone [28]." 

b) consist of negative-energy stable (i.e., stationary) bound states [29]. 

c) propagate in the Fermi sea with an energy fi^K"^ /2{2m), where m is the effective electron (or hole) 
mass and hK the CP center-of-mass momentum (CMM) [30], hence assertion # 3 above. 
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d) with the hnear dispersion E <x vphK, where is the Fermi velocity, simply represent the sound 
mode of the ideal Fermi gas (+ interactions), with sound speed vp/\fd in any dimensionality d [31]. 

e) "...with K represent states with net current flow" [32]. 

f) and BCS pairs are the same thing [33]. 

g) are not bosons, Ref. [34] p. 38, hence assertion # 2 above. 

3 Ordinary Cooper pairing 

For bosons with excitation (e.g., kinetic) energy for small CMM K given by 

SK = CsK' + o{K% (1) 

with Cs some coefficient and s > 0, BEC occurs in a box of length L if and only if d > s, since Tc = for 

all d < s. The commonest example is s = 2 as in the textbook case of ordinary bosons with exactly £k = 
fi^K^ /2mB, with ms the boson mass, giving the familiar result that BEC is not allowed for d < 2. The general 
theorem for any s > is stated as follows. The total boson number is 

Ar = Aro(r)+ ^[exp/3(£K-MB)-l]"' (2) 

with f3 = UbT. Since No{Tc) ~ while the boson chemical potential fx^ also vanishes at T = Tc, in the 
thermodynamic limit the (finite) boson number density becomes 

/■oo 

N/L'^c^Ad dK K'^-^e^^ 13 ciCsK' + ■■■)- l]-^ (3) 
where is a finite coefficient. Thus 

/•^max poo 

N/L'^ ~ Aa{kBTc/Cs) / dKK'^-'-^ +/..., (4) 

where i^max IS Small and can be picked arbitrarily so long as the integral • • • is finite, as is N/ L'^. However, 

if d = s the first integral gives hiK |^"a-=== — oo; and if d < s it gives l/(d — s)K'^^'^ |^ma,x_ _qq Hence, Tc 
must vanish if and only if d < s, but is otherwise finite. This conclusion hinges only on the leading term of the 
boson dispersion relation ek- 

The case s = 1 emerges in both the "ordinary" CP problem [29] to be recalled now, as well as in the 
"generalized" case of the next section. The question of whether or not CPs are bosons or not is resolved in 
Appendix A. For s — 1 BEC occurs for all d > 1. Striking experimental confirmation of how superconductivity 
is "extinguished" as dimensionality d is diminished towards unity has been reported by Tinkham and co- 
workers [35] [36]. They report conductance vs. diameter curves in superconducting nanowires consisting of 
carbon nanotubes sputtered with amorphous Mo'jgGe2i {Tc — 5.5 K) and of widths from 22 to 10 nm. The 
conductance diminishes to zero, implying the nanotube becomes an insulator below a certain diameter, thus 
exposing how Tc vanishes for the thinnest nanotubes. 

The CP equation for the energy £k of two fermions above the Fermi surface with momentum wavevectors 
ki and k2 (and arbitrary CMM wavenumber K where K = ki + k2) is given by 

[h^k''/m - 2Ef -£k + Tt^KyAmj^Pi, = - ^ Vkq^-q, (5) 

q 
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where k = i(ki — is the CP relative momentum and tjj^ its wave function in momentum space. The 

prime on the summation implies restriction to states above the Fermi surface with energy Ep = h'^kp/2m, viz., 
|k±K/2| > kp, and Vkq is the double Fourier transform of the interaction defined as 

Vkq = y dr J dr'e-^'i-'-y(r,r')e*--', (6) 

with V{r,r') the (possibly nonlocal) interaction in real d-dimensional space. 



3.1 Delta interaction 

If the interfcrmion interaction V{r,r') is local, then V{r,r') = V{r)S{r — r') in (6). Moreover, if y(r) = —VoS{t) 
with vo > 0, this gives Vkq = vo/L'^ and (5) becomes, for any d, 



fi2/!c2 _ 2Ef -£k + h^K"^ l^m vo 



P) 



In d = 1 where the 5- well supports a single bound state, the problem is quite tractable [37]. In either 2D or 3D, 
however, the 5- well supports an infinite set of bound levels with the lowest level in each case being infinitely 
bound. This in turn leads to a rigorous collapse of the many-fermion system [38]. To prevent this unphysical 
collapse the 5- wells must be "regularized," i.e., constructed, say, from square wells [39] such that the remaining 
(5-wcll possesses only one bound level. This leaves an infinitesimally small strength parameter vo which would 
make the rhs of (7) diverge (so as to cancel the Ihs that also diverges in 2D and 3D but not in ID). Combining 
(7) for 2D with the vacuum two-body momentum-space Schrodinger equation for the same (5-potential well 
allows eliminating [40] Vq in favor of the (positive) binding energy B2 of the single bound level of the regularized 
^-well. One arrives at 

^ B2 + fi^k'^/m " ^ /i^fcVm - 2Ep - Sk + h^K'^/4m' 

where B2 > now serves as a coupling constant. A small- ii" power-series expansion for £k gives the analytic 
expression valid for any dimensionless coupling B2/EP > 0, 



ek = £k — £0 = —hvpK - 

TT 



2\ 



l-(2-[4/7r]2) 

^2 



where a nonnegative CP excitation energy ek has been defined, and the Fermi velocity vp comes from Ep/kp = 
Hvp/2. The leading term in (9) is linear in K, followed by a quadratic term. It is clear that the leading term 
in (9) is quadratic, namely 

Ek = h^K^/2{2m) + 0{K^), (10) 

provided vp and hence Ep vanish, i.e., there is no Fermi sea. This is just the familiar nonrelativistic kinetic 
energy in vacuum of the composite (so-called "local") pair of mass 2m and CMM K. The same result (10) 
is also found to hold in 3D but not analytically as here in 2D. Figure 1 shows exact numerical results of a 
dimensionless CP excitation energy ex /{—£q) (in figure, Aq means present — f > 0) as function of K/kp for 
different couplings i?2. Note that ordinary CPs break up whenever £k turns from negative to positive, i.e., 
when £k vanishes, or by (9) when Sk/{—£o) = £k/^q = 1- These points are marked in the figure by dots. In 
addition to the exact results (full curves) also shown are some results for the linear approximation [first term 
on the right-hand side of (9), dot-dashed lines (virtually coinciding with the exact curve for all B2/EP < 0.1)], 
as well as for the quadratic approximation (dashed parabolas) as given by the leading term in (10) for stronger 
couplings. For weak enough coupling the exact dispersion relation is virtually linear — in spite of the divergence 
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of the quadratic term in (9) as B^jEp 0. As Ep decreases the quadratic dispersion relation (10) very slowly 
begins to dominate. A rcsnlt unique to 2D (and associated with the fact that in 2D the fermionic density of 

states is independent of energy) is that — £o = Ag = 

In 3D [41] instead of (8) similar procedures for two-component fermions give 

- = (11) 

k ' k,(|k±K/2l>feF) / ft -f / 

where a is the s-wave scattering length associated with the regularized (5-well, which corresponds with weak to 
strong coupling according as — oo < l/kpa < +oo. One finds for weak coupling (kpa 0~, e.g., prior to the 
well-known first-bound-state singularity as the depth of a 3D potential well is increased) that 

£o/Ef -(8/e2) exp(-7r/A;F|a|), (12) 

a result first reported by Van Hove [42]. For strong coupling {kpa 0+, beyond the single-bound-state 
resonance) one gets 

£o/Ef ^ ~2/{kFaf. (13) 

Numerical results in 3D very similar to those in Fig. 1 for 2D are obtained. Namely, for weak coupling the CP 
dispersion curves are very nearly linear while for smaller density they very slowly tend to the quadratic. The 
limit given by (9) in 2D was found too complicated in 3D to be evaluated analytically, except for weak coupling. 
Repeating the 2D analysis without attempting to explicitly determine the coefficient of the quadratic term one 
gets, for the (5-well interfermion interaction in weak coupling, 

eK = £K-£o = \hvFK + 0{K^). (14) 

This is the same result cited without proof by Schrieffer in 1964 (Ref. [34] , p. 33) for the BCS model interaction, 
to which we now turn. 



3.2 BCS model interaction 

The BCS model interaction is of the simple form 

{ -V if fcf <|k±iK|, |k'±iK|< + 
I otherwise, 

where Vkk' is defined in (6). Here F > 0, and hu)D = h'^kjj/2m is the maximum energy of a vibrating-ionic- 
lattice phonon. This means that two fermions interact with a constant attraction —V when the tip of their 

relative-momentum wavevector k points anywhere inside the overlap volume in fc-space of the two spherical 
shells in Figure 2. Inserting (15) into (5) and converting sums over k into energy integrals by introducing the 
electronic density of states (DOS) g{e) gives 

^From this one immediately obtains the familiar result for K = 0, 

Here A = g{EF)V is a dimensionless coupling constant and giEp) the DOS for each spin evaluated at the Fermi 
energy. The equality in (17) is exact in 2D for all coupling — as well as in ID or 3D provided only that huD 
<C Ef so that g{€) ~ g{EF), a constant that can be taken outside the integral in (16). 
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For a 2D system (16) gives [43] for weak coupling 

£k So + (2/TT)hvFK + OiK'^). (18) 

The exact dispersion relation obtained numerically from (16) for A = ^ and Hlod/Ep = 10~^ shows that the 
linear approximation (18) is very good for moderately small A and tuvo/Ep, over the entire range of K values 
for which £k < 0. Note that the linear term carries the same coefficient as (9) for a different interfermion 
interaction. Pair breakup, specifically £k > for these values of A and fuvD/Ep, occurs at a relatively small 
value of K, about four orders of magnitude smaller than the maximum value + k'jj allowed by the 

interaction (15). 

In 3D, assuming Huid/Ep << 1 so that the 3D DOS g{e) = {L'^/n'^h'^)y/ni^e/2 can be replaced by g(Ep) 
and then taken outside the integral sign, the result cited in Ref. [34], p. 33, (see also Ref. [44], p. 336, Prob. 
10.4 but note here an erroneous coeffcient) follows, namely 

Sk So + ihvpK + OiK^). (19) 

Exact numerical results in 3D are qualitatively similar [43] to those in 2D as regards goodness of the linear 
approximation for weak coupling. 



4 Generalized Cooper pairing 

The "ordinary" CP problem just summarized for two distinct interfermion interactions (the (5- well and the BCS 
model interaction) neglects the effect of two-hole (2h) CPs treated on an equal footing with two-particle (2p) 
CPs — as Green's functions [44] can naturally guarantee. On the other hand, the BCS condensate consists [13] 
of equal numbers of 2p and 2h Cooper " correlations This was already evident, though scarcely emphasized, 
from the perfect symmetry about /x, the electron chemical potential, of the well-known Bogoliubov [45] u^(e) 
and it^(e) coefficients [see just below (30) later on], where e is the electron energy, the prime motivation here 
rests on the recently established remarkable fact [13] that a BCS condensate is precisely a BE condensate with 
equal numbers of 2p and 2h zero-CMM CPs, in the limit of weak coupling. Further motivation comes from 
the unique but unexplained role played by hole charge carriers in the normal state of superconductors [46] in 
general (see also Ref. [47]). Final motivation stems from the ability of the "complete (in that both 2h- and 
2p-CPs are allowed in varying proportions) BF model" of Refs. [13]- [15] to "unify" both BCS and BEC theories 
as special cases, and to predict substantially higher Tc's than BCS theory without abandoning electron-phonon 
dynamics. The latter is important as compelling evidence for a significant, if not sole, presence of it in high-T^ 
cuprate superconductors from angle- resolved photoemission spectroscopy data has recently been reported [48] . 

In this section we sketch how the Bethc-Salpeter (BS) many-body equation (in the ladder approximation) 
treating both 2p and 2h pairs on an equal footing shows that, while the ordinary CP problem [based on an ideal 
Fermi gas (IFG) ground state (the usual "Fermi sea")] does not possess stable energy solutions, but that it does 
so when the IFG ground state is replaced by the BCS one. This is equivalent to starting from an unperturbed 
Hamiltonian that is the BCS ground state instead of the pure-kinetic-energy operator corresponding to the 
IFG. We discuss how: i) CPs based not on the IFG-sea but on the BCS ground state survive in a nontrivial 
solution as "generalized" or "moving" CPs which are positive energy resonances with an imaginary energy term 
leading to finite-lifetime effects; ii) as in the "ordinary" CP problem of the previous section, their dispersion 
relation in leading order in the total (or center-of-mass) momentum (CMM) UK = /l(k^ -|- k2) is also linear 
rather than the quadratic h'^K'^ /2{2m) of a composite boson (e.g., a deuteron) of mass 2m moving not in the 
Fermi sea but in vacuum; and iii) this latter "moving CP" solution, though often confused [31] with it, is 
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physically distinct from another more common trivial solution sometimes called, even though Bogoliubov [8] 

was the first to derive it, the Anderson-Bogoliubov-Higgs (ABH) [49], ([9] p. 44), [50]-[52] collective excitation. 
The ABH mode is also linear in leading order and goes over into the IFG ordinary sound mode in zero coupling. 
All this occurs in both 2D [53] as well as in the 3D study outlined earlier in Ref. [54]. We focus here on 
2D because of its interest [55]-[57] for quasi-2D high- Tc cuprate superconductors. In general, the results will 
be crucial for BEG scenarios employing BF models of superconductivity, not only in exactly 2D as with the 
Bcrczinskii-Kostcrlitz-Thoulcss [58] [59] transition, but also down to (1 + e)D which characterize the quasi-lD 
organo-metallic (Bechgaard salt) superconductors [60]- [62]. These results also apply, albeit with a different 
interaction, to neutral-atom superfluidity as in liquid "^He [63] and very probably also in ultracold trapped alkali 
Fermi gases such as [64] and ^Li [65] atoms. 

In dealing with the many-electron system we again assume the BCS model interaction (6) in 2D with double 
Fourier transform 

i'{ki,k[) = -{k%/kik[)V ii kp - kn < ki,k[ < kp + kn, 

= otherwise. (20) 

Here y > 0, hkp = mvp the Fermi momentum, m the effective electron mass, vp the Fermi velocity, and 
ko = ujd/vf with w^) the Debye frequency; note difference with previous definition just below (15). The usual 
physical constraint huJo ^ Ep then implies that kp,/kp = fuvD/^Ep <C 1. Assuming perfect ph symmetry 
about the Fermi surface, we set 

ekC^Ep + hvp{k - kp) (21) 
as it simplifies all calculations when very near the Fermi surface. 

The bound-state BS wavefunction equation [54] in the ladder approximation with both particles and holes 
for the original IFG-based CP problem is 

*(k,£;) = -{if Go (K/2 + k, £kI2 + E) Go (K/2 - k, £k/2 - E) x 

/ di; J^Ek'^(|k-k'|)*(k',E')- (22) 

— oo 

Here L'' is the "volume" of the d-dimensional system; K = ki -|-k2 is the CMM and k = ^(k^ — k2) the relative 
wavevectors of the 2e bound state whose wavefunction is ^{k,E); £k = Ei + E2 is the energy of this bound 
state while E = Ei — E2, and Go {1^/2 + \c,£/2 + E) is the bare one-fermion Green's function given by 

0^1' ^' i\-Ei + eu,-Ep-ie -Ei + eu, - Ep + is } ^ ' 

where eki = ti?k\/2m and ^(a;) = 1 for a; > and = for a; < 0, so that the first term refers to electrons 
and the second to holes. Figure 3 shows all Feynman diagrams for the 2p. 2h and ph wavefunctions 
'0_ and respectively, that emerge in the general (BGS-gromid-state-based) problem to be discussed later. 
For the present IFG-based case, diagrams in shaded rectangles do not contribute as they involve factors of 
6{k\ — kp)d{kp — ki) = 0. Since the energy dependence of ^'(k, E) in (22) is only through the Green's functions, 
the ensuing energy integrals (22) can be evaluated directly in the complex £'-plane and yield, for interaction 
(20), an equation for the wavefunction ^j,. in momentum space for GPs with zero GMM K = ki + k2 = that 
is 

(2e,-£:o)V'k = V^E'^k'-V^E"^k'- (24) 

k' k' 

Here = h'^k'^/2m — Ep while £0 is the K = eigenvalue energy, and k = |(ki — k2) = k^. The single 
prime over the first (2p-GP) summation term denotes the restriction < < hcun while the double prime 
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in the last (2h-CP) term means —Hcod < ^k' < 0- Without this latter term we have Cooper's Schrodinger-like 
equation [29] for 2p-CPs whose implicit solution is clearly ipj^ = (2^^, — £a)~^V V'k'- Since the summation 
term is constant, performing that summation on both sides allows canceling the ^/;j(.-dependent terms, leaving 
the eigenvalue equation X)ij(2^j. — £o)~^ = l/V. This is one equation in one unknown Sq; transforming the 
sum to an integral over energies gives the familiar solution (17) exact in 2D, and to a very good approximation 
otherwise if hu)D <C Ep, where A = VN{Ef) with N{Ef) the electronic DOS for one spin. This corresponds 
to a negative-energy, stationary-state bound pair. For K the 2p-CP eigenvalue equation becomes 

5^'(2Cfc -£k + n^K^/Am)-^ = 1/V. (25) 

k 

Note that a 2p CP state of energy £x is characterized only by a definite K but not definite k, in contrast to a 
"BCS pair" defined [Ref. [7], Eqs. (2.11) to (2.13)] with fixed K and k (or equivalently definite ki and k2); see 
Appendix A. Without the first summation term in (24) the same expression (17) for Sq for 2p-CPs follows for 
2h-CPs, apart from an overall sign change. 

The complete equation (24) cannot be derived from an ordinary (non-BS) Schrodingcr-like equation in spite 
of its simple appearance. To solve it for the unknown energy £q, let the rhs of (24) be defined as A — B, with 
A relating to the 2p term and B to the 2h one. Solving for the unknown tp^ gives 

V'k = (^ - S)/(2Cfc - So) or equivalently V(0 = - B)/{2^ - So) (26) 

whence 

A = \ d£,ip{£,)^ -{A- B)\ dz/z={A-B)x (27) 

Jo 2 J_£g 

B = \ f d^^P{0 = l{A-B)X [ dz/z={A-B)y. (28) 

The integrals are readily evaluated giving x = iAln(l — 2hiUD/So) and y = — iAln(l + 2fujjD/So)- As A and B 
still contain the unknown ^p{^) let us eliminate them. Note that (27) and (28) are equivalent to two equations 
in two unknowns A and B, namely 

{l-x)A + xB = 
-yA+{l+y)B = 0. 

This leads immediately to the equation 1 — x + y — 0, which on inserting the definitions of x and y becomes 

l = ^X\n[l-{2nojD/So)% 

which finally yields 



So = ±i2nu;D/Vey^ - 1. (29) 

As the CP energy is pure- imaginary there is an obvious instability of the CP problem when both particle- and 
hole-pairs are included. This was reported in Refs. [9] p. 44 and [66], who did not, however, stress the pure 
2p and 2h special cases just discussed. Clearly then, the original CP picture is meaningless if particle- and 
hole-pairs are treated on an equal footing, as consistency demands. 



5 BCS-based BS treatment of Cooper pairing 

However, a BS treatment not about the IFG sea but about the BCS ground state vindicates the CP concept as 
a nontrivial solution. This is equivalent to starting not from the IFG unperturbed Hamiltonian but from the 
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BCS one. Its physical justification lies in recovering three expected items: ABH sound mode, the BCS T = 

gap equation and finite-lifetime efi^ccts of the "moving CPs." In cither 2D [53] or 3D [54] the BS equation yields 
a 4 X 4 determinant which reduces to a 3 x 3 and a 1 x 1 determinant representing, respectively, two distinct 
solutions: a) the trivial ABH sound solution and b) a highly nontrivial moving CP solution, respectively. In 
either case the BS formalism gives a set of three coupled equations, one for each (2p, 2h and ph) channel 
wavefunction for any spin-independent interaction such as (20). However, the ph channel decouples, leaving 
only two coupled wavefunction equations for the ABH solution in 2D which we consider first. We note that in 
Ref. [67] the hh channel was ignored, leading to a 3 x 3 determinant from which only the trivial ABH solution 
emerges, but the nontrivial moving CP one was missed entirely. Thus, the IFG Green function (23) is replaced 
by the BCS one 

where = -y/^^TTA^ with A the T = fermionic gap, vl = ^{l — ^k/^k) and ul = 1 — are the 
Bogoliubov functions [68]. As A ^ these three quantities become 9{ki — kp) and 0{kp — fci), respectively, 
so that (30) reduces to (23), as expected. Substituting Go(ki,i?i) by Go(ki,i?i) corresponds to rewriting the 
total Hamiltonian so that the pure-kinetic-energy unperturbed Hamiltonian is replaced by the BCS one. The 
remaining Hamiltonian terms are then assumed suitable to a perturbation treatment. We focus in this section 
only on 2D. 



5.1 ABH sound (trivial) solution 



The equations involved are too lengthy even in 2D and will be derived in detail elsewhere, but for the trivial 
ABH sound solution the aforementioned 3x3 determinant boils down to the single expression 



— XhvF / dk #{WK/2+k'"K/2-k + 'l^K/2+k1^K/2-k}x 

^'^ JkF-ko Jo 



VK/2+kVK/2-k ^K/2+k ^K/2-k ^ _ ^ , . 



£k + + -£^K/2-k -^Jf + -E'K/2+k + -2^K/2-k" 

where <^ is the angle between K and k. Here ko = ujd/vf', note difference with definition just below (15). 
As before A = VN{Ep) with N{Ef) = the constant 2D electronic DOS and F > is the interaction 

strength defined in (20). Angle- resolved photocmission spectral studies of BiSrCaCuO have shown evidence 
[69] in this cuprate for the Bogoliubov functions u'^ and f ^ , both above and below the Fermi energy. 

The ABH collective excitation mode energy £k must then be extracted from this equation. For K = it is 
just So = (Ref. [9] p. 39). Then (31) rewritten as an integral over ^ = h^k'^/2m — Ep reduces to 

rhwo I 

/ d^/^^e + A2 = 1/A, (32) 

JO 

or the familiar BCS T = gap equation for interaction (20). The integral is exact and gives 

A= ;zwr,/sinh(l/A). (33) 
Returning to the ABH energy £k equation (31) and Taylor-expanding £k about K = Q and small A leaves 

£k = ^K + 0{K^) + o{\), (34) 
v2 
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where o(A) denote interfermion interaction correction terms that vanish as A ^ 0. Note that the leading term 
is just the ordinary sound mode in an IFG whose sound speed c = Vp/Vd in d dimensions. 

The latter also follows elementarily on solving for c in the familiar thermodynamic relation dP/dn = m(? 
involving the zero-temperature IFG pressure 

P = n'^[d{E/N)ldn] = 2nEF/{d+2) 

= 2Cdn2/<*+V(d + 2) (35) 

where the constant Ca will drop out. Here the IFG ground-state (internal) energy per fermion E/N = dEp/{d+ 
2) = Cdn'^/'^ was used along with Ep = h'^k%/2m and 

n = N/L'^ = kj, /d2'^-\^/'^T{d/2) (36) 

for the fcrmion-mimbcr density n. The derivative of (35) with respect to n finally gives c = Hkp /rnVd = vp/Vd, 
which in 2D is just the leading term in (34). 



5.2 Moving CP (nontrivial) solution 

The nontrivial moving CP solution of the BCS-ground-state-based BS treatment, which is entirely new, comes 
from the remaining 1x1 determinant. It leads to the pair energy 6k which in 2D is contained in the equation 



— XhVp / dk #WK/2+k^^K/2-kX 

JkF-ko Jo 



27r 

X |UK/2-k^'K/2+k - MK/2+k^^K/2-k| c2 , ^ f N2 = (-^T) 

-t^K + (■C'K/2+k + -C'K/2-kj 

In addition to the pp and hh wavefunctions (depicted diagrammatically in Ref. [54] Fig. 2), diagrams associated 
with the ph channel give zero contribution at T = 0. A third equation for the ph wavefunction describes the 
ph bound state but turns out to depend only on the pp and hh wavefunctions. Taylor-expanding £k in (37) in 
powers of K around K = 0, and introducing a possible damping factor by adding an imaginary term —iVx in 
the denominator, yields to order K"^ 

^ , 1 ^Vp ^i/A^^2 



±£k o,2A+ —hvpK + -^e^/^K' 
2iT 9 ko 



— I 



^nvpK+^^e'/'K' 
TT 12 ko 



+ 0{K^) (38) 



where the upper and lower signs refer to 2p- and 2h-CPs, respectively. A linear dispersion in leading order 
again appears, but now associated with the bosonic moving CP. Hence the positive-energy 2p-CP resonance has 
a width Vk and a lifetime 

TK = h/2TK = h/2 \{\lTJ:)hvpK + {hvp/l2kD)e^'^K'^\ ■ (39) 

This diverges only at if = 0, falling to zero as K increases. Thus, "faster" moving CPs are shorter-lived 
and eventually break up, while "non-moving" ones are infinite-lifetime stationary states. The linear term 
{X/2TT)hvpK in (38) contrasts sharply with the coupling-independent leading-term in (18) [or 1/2 in 3D in (19), 
Ref. [34] p. 33, instead of 2/-jt] that follows from the original CP problem (25) neglecting holes. This holds 
for either interaction (20) [43] or for an attractive delta interfermion potential well in 2D [40] or in 3D [41]. 
These 5-wells are imagined regularized [39] to possess a single bound level whose; binding energy (in 2D) or 
scattering length (in 3D) serve as the coupling parameter. Figure 4a shows the exact moving CP (mCP) energy 
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(full curves) extracted from (37), along with its leading linear-dispersion term (thin short-dashed lines) and this 

plus the next (quadratic) term (long-dashed curves) from (38). The interaction parameter values used with (20) 
were Hlod/Ef = 0.05 (a typical value for cuprates) and the two values A = j and ^. Using (33) in (38) gives 

£o/Ef = 2A/Ef = 2hiJD/EFsmh{l/X), (40) 

having the values ~ 0.004 and 0.028, respectively (marked as dots on the figure ordinate). Remarkably enough, 
the linear approximation (thin short-dashed lines in figure) is better over a wider range of K/kp values for weaker 
coupling (lower set of three curves) in spite of a larger and larger (because of the factor e^^^) partial contribution 
from the quadratic term in (38). This peculiarity also emerged from the ordinary CP treatment of Sec. 3, Refs. 
[40] [41] [43] . It suggests the expansion in powers of K to be an asymptotic series that should be truncated after 
the linear term. For reference we also plot the linear leading term hvpK/ \f2 of the sound solution (34). We note 
that the coupling-independent leading-term [43] {2/'K)hvFK from the original CP problem neglecting holes, if 
graphed in Fig. 4, would almost coincide with the ABH term hvFK/^/2 and have a slope about 90% smaller. 

Empirical evidence for the linearly- dispersive nature of Cooper pairs in cuprates has been argued by Wilson 
[70] to be suggested by the scanning tunneling microscope conductance scattering data in BSCCO obtained by 
Davis and coworkers [71] [72]. More suggestive direct evidence is shown in Figure 5 [73] with experimental data 
(mostly from penetration-depth measurements) for two 3D SCs [74] [75], two quasi-2D cuprates [76]-[78], and 
a quasi-lD SC [79]. The data are seen to agree quite well, at least for T > 0.5Tc, with the pure-phase (only 
2e- or 2h-CP) BEC condensate fraction formula 1 — (T/TcY^'' for d = 3, 2 and 1, respectively, provided one 
assumes s = 1. For lower T's, one can argue on the basis of Fig. 9 below that a mixed BEC phase containing 
both 2e- and 2h-CPs becomes more stable (i.e., has lower Helmholtz free energy) so that the simple pure-phase 
1 — (T/TcY^^ formula is no longer strictly valid. 

As in Cooper's [29] original equation (25), the BS moving CPs are characterized by a definite K and not 
also by definite k as the pairs discussed by BCS [7]. Hence, the objection does not apply that CPs are not 
bosons because BCS pairs with definite K and k (or equivalently definite ki and k2) have creation/annihilation 
operators that do not obey the usual Bose commutation relations [Ref. [7], Eqs. (2.11) to (2.13)]. In fact, 
either (25) or (37) shows that a given "ordinary" or "generalized" CP state labeled by either K or £k can 
accommodate (in the thermodynamic limit) an indefinitely many possible BCS pairs with diflEerent k's; see 
Ref. [82] or Appendix A. A recent electronic analog [83] of the Hanbury Brown- Twiss photon-effect experiment 
suggest electron pairs in a SC to be definitely bosons. 

To conclude this section, hole pairs treated on a par with electron pairs play a vital role in determining the 
precise nature of CPs even at zero temperature — only when based not on the usual IFG "sea" but on the BCS 
ground state. Their treatment with a Bethe-Salpeter equation gives purely-imaginary-energy CPs when based 
on the IFG, and positive-energy resonant-state CPs with a finite lifetime for nonzero CMM when based on the 
BCS ground state — instead of the more familiar negative-energy stationary states of the original IFG-based CP 
problem that neglects holes, as sketched just below (24). The BS "moving-CP" dispersion relation (38), on the 
other hand, resembles the plasmon dispersion curve in 3D. It is gapped by twice the BCS energy gap, followed 
by a linear leading term in the CMM expansion about K = 0, instead of the quadratic for the 3D plasmon 
curve. This linearity is distinct from the better-known one (34) associated with the sound or ABH collective 
excitation mode whose energy vanishes at K = 0. Thus, BF models assuming this CP linearity for the boson 
component, instead of the quadratic h^K"^ /2{2m) assumed in Refs. [3], [11]-[14], [89]-[92] among many others, 
can give BEC for all d > 1, including exactly 2D. Such BF models can then in principle address not only 
quasi-2D cuprate but also quasi-lD organo-metallic superconductors. 
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6 The CBFM Hamiltonian 



The CBFM [13, 14] is described (in d dimensions) by the Hamiltonian H = Hq + Hint- The unperturbed 
Hamiltonian Hq corresponds to a "normal" state which is an ideal (i.e., noninteracting) gas mixture of unpaired 
fermions and both types of CPs, two-electron (2e) and two-hole (2h), and is given by 

Ho=Yl ek.a+ _,^ak„., +J2E+{K)b+bj^ -Y,E_{K)c+ck (41) 

ki,si K K 

where as before K = ki + k2 is the CP CMM wavevector, k = ^(kj^ — k2) its relative wavevector, while 
Cki = h?kf/2m are the single-electron, and E±{K) the 2e-/2h-CP phenomenological, energies. Here a^i si 
(<iki,si) are creation (annihilation) operators for fermions and similarly 6^ (6k) and (ck) for 2e- and 2h-CP 
bosons, respectively. Two-hole CPs are considered distinct and kinematically independent from 2c-CPs since 
their Bose commutation relations involve a relative sign change, in sharp contrast with electron or hole fermions 
whose Fermi anticommutation relations do not. Our present formulation is of course nonrelativistic because of 
relatively low temperatures. 

At the opposite extreme of very high T's (compared with the boson rest-mass energy), however, one has 
a relativistic regime where pair production takes place and BEC must take N antibosons of charge, say, —q 
into account along with the N bosons of charge q. In units such that h = c = ks = I the boson energy 
is ek = {K^ + m\Y^^. Charge conservation requires that not N = No{T) + X^k^^o [®-^P /^(^^ ~ I-'-b) ~ 
be constant but rather N— N, where N is the same expression as A'' but with instead of —^b- If 

p = q{N — N)/L^ = qn is the net conserved charge density, it is shown in Ref. [93] that Tc = (3]n]/ms)5 and 
that the condensate fraction na{T)/n = 1-{T/Tcf. This is qualitatively different from the better-known results 
assuming only N constant, which arc the mass-independent Tc = [7r2n/C(3)]^/^ and no{T)/n = 1 - {T/Tcf. 
This analogy with the CBFM exhibits the strikingly dramatic effect of including or not antiparticles — or, in our 
nonrelativistic problem, of including or not hole pairs. 

The interaction Hamiltonian Hint consists of four distinct BF interaction vertices, see Fig. 6, each with 
two-fermion/one-boson creation or annihilation operators, depicting how unpaired electrons (subindex -|-) [or 
holes (subindex — )] combine to form the 2e- (and 2h-CPs) assumed in the d-dimensional system of size L, 
namely 

Hint=L-'^'^^f+{k){a+^^^^^a+^^,^^^h^ + a-k-FiK,i«k+iK,T^K} 

k,K ^ ' ^ ' 

+i"''^'IZ/-(fc){ak+iK,T"-k+iK,i^K + a_k+iK,i«k+iK,TCK}- (42) 

k,K ^ ' ^ ' 

Note that the fermion-pair interaction Hint is reminiscent of the Frohlich (or Dirac QED) interaction Hamilto- 
nian (Ref. [44], p. 396 ff.) involving two fermion and one boson operators, but with two types of CPs instead 
of phonons (or photons). But in contrast with Frohlich and Dirac there is no a conservation law for the number 
of unpaired electrons, i.e., [Hint, ^kiflki.s "^ki.sj 7^ 0. (Note however that [Hint, ^ ^^i^ki.s ^ki.sj = 

ki,si ki,si 

and [Hint, 5^ *i'^ki , aki,sj — 0.) Just as the Frohlich (or Dirac) interaction Hamiltonians are the most 

ki,.si ' ^ 

natural ones to use in a many-electron/phonon (or photon) system, one can conjecture the same of (42) for 
the BF system under study. Indeed, this Hint has formally already been employed under various guises by 
several authors [10]-[12],[84]-[87], but without hole pairs. More recently, a similar BF Hint has been employed 
[89]- [92] to study quantum degenerate Fermi gases consisting of neutral '^'^K atoms and their so-called Feshbach 
"resonance superfluidity." However, these authors assume quadratically-dispersive CPs, besides also excluding 
2h-CPs and cannot thus relate their formalism to BCS theory. 
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The energy form factors f± (fc) in (42) are essentially the Fourier transforms of the 2e- and 2h-CP intrinsic 
wavefunctions, respectively, in the relative coordinate between the paired fermions of the CP. In Refs. [13] [14] 
they are taken simply as 

/.(.) = { i '»'^' <;<^/+*=. (43) 
I otherwise, 

/-(e) = ii f-^r-*^<-%. ,44) 
I otherwise. 

The quantities Ef and 5e are new phenomenological dynamical energy parameters (in addition to the positive 
BF vertex coupling parameter /) that replace the previous such E±{0), through the relations 

Ef = i [E+ (0) + E_{0)] and 5s = ^[E+{0) - E_{0)], (45) 

where E±{0) are the (empirically unknown) zero-CMM energies of the 2e- and 2h-CPs, respectively. Clearly 

E±{0) = 2Ef ± Sc. The quantity Ef will serve as a convenient energy scale and is not to be confused with the 
Fermi energy Ep — ^mvp = ksTp where Tp is the Fermi temperature. The Fermi energy Ep equals nh^n/m 
in 2D and (fi^/2m)(37r^n)^/'^ in 3D, with n the total number-density of charge-carrier electrons. The quantities 
Ef and Ep coincide only when perfect 2e/2h-CP symmetry holds. 

The interaction Hamiltonian (42) can be further simplified by keeping only the K = terms. One can then 
apply the Bogoliubov "recipe" [68] (see, also e.g., [44] p. 199) of replacing in the full Hamiltonian H = Ho + Hint 
all zero-CMM 2e- and 2h-CP boson creation and annihilation operators by their respective c-numbers ^/N^ and 
y/Mo, where No{T) and Mo{T) are the number of zero-CMM 2e- and 2h-CPs, respectively. One eventually seeks 
the lowest critical temperature Tg such that, e.g., No{Tc) or Mo(Tc) vanish. Note that Tc calculated thusly can, 
in principle, turn out to be zero, in which case there is no EEC. Now, H — /iTV can be diagonalized exactly via 
a Bogoliubov- Valatin transformation [45] in terms of the thermodynamic potential f2 = —PL'^ for the CBFM, 
with P its pressure and L'' the system "volume," which is 

Q(T, L'^, /i, Nq, Mo) = -keTln rTre"'''^-''*)] , (46) 

where "Tr" stands for "trace." Inserting (41) and (42) into (46) one obtains [13] after some algebra 

n{T, /i, iVo, Mo)/L'' = 

= deN{e)[e- fi- E{e)]-2kBT deN{e) ln{l + exp[-P E{e)]}+ 
Jo Jo 

+ [E+{0) - 2/i] no + keT / drjMir]) ln{l - exp[-/?f+(r?)]}+ 

Jo+ 

POO 

+ [2n- E_(0)] mo + keT / drjMir]) ln{l - exp[-/Jf_(ry)]}. (47) 

Jo+ 



For = 3 one has 



m3/2 _ 2m3/2 



for the (one-spin) fermion and boson DOS at energies e = h?k'^/2m and jj = h^K'^ /2{2m), respectively. The 
latter is an assumption to be lifted later so as to include Fermi sea effects which transform the boson dispersion 
relation from quadratic to linear, as already discussed in Sections 3 to 5. Also, recalling the 2e- and 2h-boson 
energies E± (K) introduced in (41) the (bosonic) energies £+ (jy) and £- {rj) to simplify notation have been defined 
as 

E+{K) -2^I = E+{0)-2^i + 7J = £+{r,), (49) 
2ii- E_{K) =2ii- E_{0) +Tj= £_{tj). (50) 
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Finally, the relation between the fermion spectrum £'(e) and fermion energy gap A(e) is of the form 

E{e)= v^(e-^)2 + A2(e), (51) 
A(e) = V^/+(e) + V^/_(e). (52) 

This last expression for the gap A(e) implies a simple r-clepcndence rooted in the 2e-CP no(T) = Nq{T)/L'^ 
and 2h-CP mo{T) = Mq{T)/L'^ number densities of BE-condensed bosons, i.e., 

A(T) = v/;^T)/+(e) + v/;;^T)/_(e). (53) 

A A(T)^ temperature dependence in the order parameter in thin films of the Tl2Ba2CaCu20e. cuprate suggests 
itself from thermal difference reflectance spectra [95]. This implies that condensate densities no(T) or mo{T) 
might be more fundamental as order parameters than A(T), at least for this material. 

Minimizing il{T, L"^, fi, Nq, Mq) with respect to Nq and Mo, and simultaneously fixing the total number N 
of electrons by introducing the electron chemical potential /z, one thus specifies an equilibrium state of system 
with volume L'^ and temperature T by requiring that 

0, =0, and = -N. (54) 



dNo ' dMo ' diJ, 

Here N evidently includes both paired and unpaired CP fermions. Some algebra then leads to the three coupled 
transcendental Eqs. (7)-(9) of Ref. [13]. These can be rewritten somewhat more transparently as: a) two 
"gap-like equations" 

oo 

2^[E+iO) - 2/.] = J de7V(e)^M^ tanh i/3E(e); (55) 



oo 

2v^[2m - E.{0)] = J deiV(e)^i|^ tanh \pE{e), (56) 



and b) a single "number equation" (that ensures charge conservation) 

2nB (T) - 2mB (T) + n/ (T) = n. (57) 

Here n = N/L^ is the number density of electrons while nB(T) and ms(T) are the number densities of 2e- and 
2h-CPs in all bosonic states, respectively, and are given by 



nBiT) ^ no(T) + / drjM{7j) ^^^^^^^ _ ^ , (58) 



0+ 



msiT) = mo{T) + j driM{r^)-^^-^^. (59) 

Clearly, 



0+ 



nf{T) ^ I deN{e)[l - ^ tanh ^f3E{e)] (60) 



is the number density of unpaired electrons at any T. Self-consistent (at worst, numerical) solution of the three 
coupled equations (55) to (57) yields the three thermodynamic variables of the CBFM 

no{T,n,ii), mo{T,n,ii), and n{T,n). (61) 
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The pressure P, entropy S and specific heat at constant volume C of an equilibrium state characterized by 
T and n are then given by 

PiT,n) = -n/L^, S{T,n)/L'' = -kB^{^lL% (62) 
C{T,n)/L'' = T-^[S{T,n)/L''], (63) 
all evaluated at fixed no{T, fj.,n), mo{T,fj.,n) and n{T,n). The Helmholtz free energy 

F{T, L'^,N) = E-TS=n + fiN (64) 
where E is the internal energy, then follows from 

F{T,n)/L'^ = -P{T,n) + niJ.{T,n). (65) 

Finally, the critical magnetic field is 

H^{T, n)/8 7r = F„(r, n)/L'' - F,(T, n)/L'' = 

= Ps (T, n) - P„ (T, n) + [m„ (T, n) - m« (T, n)] n, (66) 
with subscripts s and n meaning "superconducting" and "normal" states. 

7 Main statistical theories as special cases of CBFM 

Most significantly, the three CBFM equations just stated contain five different theories as special cases, see flow 
chart in Fig. 7. Perfect 2e/2h CP symmetry nsiT) = TOb(T) and no{T) = mo{T) can be seen from (55) and 
(56), as well as from E±{0) = 2Ef ± 5s, to imply that Ef coincides with /x. The CBFM then reduces to: 

i) the gap and number equations of the BCS-Bose crossover picture [96] for the BCS model interaction — 
if the BCS parameters V and Debye energy Hlod are identified with the BF interaction Hamiltonian Hint 
parameters f^/2Se and Se, respectively. The crossover picture for unknowns A(T) and /i(T) is now supplemented 
by the key relation 

A(r) = /v'no(T) = /x/mo(T). (67) 

The 35-year-old crossover picture is associated with names such as (in chronological order) Friedel and co- 
workers [97], Eagles [98], Leggett [99], Miyake [100], Nozieres [101], Micnas, Ranninger & Robaszkiewicz [10], 
Randeria [102], van der Marel [103], Bar- Yam [85], Drechsler & Zwerger [104], Hausmann [105], Pistolesi & 
Strinati [106], and many others. However, it seems to be a very modest improvement over BCS theory since 
the unphysically large A of about 8 is required to bring fi{Tc)/Ep down from 1.00 to 0.998 [107]. If one imposes 
that n{Tc) — Ep exactly, as follows for weak BF coupling / from the number equation, the crossover picture is 
well-known to reduce to: 

ii) ordinary BCS theory which is characterized by a single equation, the gap equation for all T. Thus, 
the BCS condensate is precisely a BE condensate whenever both nsiT) = mB{T) and no(T) = mo(T) and the 
BF coupling / is small. Indeed, for small coupling A the CBFM T = superconducting state has a lower energy 
than the corresponding BCS state, see Appendix B. The BCS state comes from a variational trial wave function 
so that its energy expectation value by the Rayleigh-Ritz variational principle is a rigorous upper bound to the 
exact value. Thus, the CBFM has a somewhat larger condensation energy than that of BCS theory, but both 
coincide in leading order. In addition, at least two universal constants, the gap-to-Tc ratio 2A(0)/fcsTc and the 
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specific heat jump ACs(Tc) to the normal value C„(Tc), have been shown [14] to coincide with the BCS values 
of 3.53 and 1.43 in this limit of the CBFM. 

On the other hand, for no 2h-CPs present the CBFM reduces [13] also to: 

iii) the BEC BF m,odel in 3D of Friedberg and Lcc [11, 12] characterized by the relation A(T) = 
/■\/no(T). With just one adjustable parameter (the ratio of perpendicular to planar boson masses) this theory 
fitted [12] cuprate T^/Tp data quite well. The ratio turned out to be 66,560 — ^just under the 10^ anisotropy 
ratio reported [108] almost contemporaneously for BSCO. Finally, when / = this model reduces to: 

iv) the ideal BF model (IBFM) of Ref. [17, 18] that predicts nonzero BEC TcS even in 2D. The "gapless" 
IBFM cannot describe the superconducting phase. But considered as a model for the normal state it should 
provide feasible TcS as singularities within a BE scenario that are approached from above Tc. Figure 8 displays 
the Tc predictions [18] in 2D for cuprate superconductors with no adjustable parameters. Finally, the CBFM 
reduces to: 

v) the familiar Tc-formula of ordinary BEC in 3D, albeit as an implicit equation with the boson number- 
density being T-dependent. 

8 Enhanced TcS from the CBFM 

The very general CBFM has been applied and gives sizeable enhancements in TcS over BCS theory that emerge 
for moderate departures from perfect 2e/2h-pair symmetry. This is attained for the same interaction model 
(with the same coupling strength A and cutoff Hujd parameters) used in conventional superconductors. The 
BF (two-fermion) interaction (42) with (43) and (44) bears a one-to-one correspondence with the more familiar 
"direct" four-fermion electron-phonon interaction, mimicked, e.g., by the BCS model interaction. Its double 
Fourier transform is a negative constant — V, nonzero only within an energy shell 2huJD about the Fermi surface, 
with uiD the Debye frequency. The correspondence is realized [13, 14] by setting /^/2fc = V and Se = Hujd, 
from which the familiar dimensionless BCS model interaction parameters A = N{Ef)V and hoJo/Ep are then 
recovered. 

The three coupled equations (55) to (57) of the CBFM that determine the rf-dimensional BE-condensate 

number-densities r7,o(T) and m.Q{T) of 2c- and 2h-CPs, respectively, as well as the fcrmion chemical potential 
/i(T), were solved numerically in 3D for A = 1/5 and huj]j/Ep = 0.001 in Ref. [14] assuming a quadratic boson 
dispersion relation rj = fi^K^ /2{2m), i.e., s = 2 in (1). For this case Figure 9 maps the phase diagram in the 
vicinity of the BCS Tc value (marked BCS-B in the figure) at An = n/n/ — 1 = (corresponding to perfect 
2e/2h-CP symmetry) where n/ is a very special case of n/(T) as defined in (60). Namely, n/ can be seen [14] to 
be the number of unpaired electrons at zero coupling and temperature. Besides the normal phase (n) consisting 
of the ideal BF gas described by Hq, three different stable BEC phases emerged plus several metastable, i.e., of 
higher Helmholtz free energy. These are two pure phases of either 2e-CP (s+) or 2h-CP (s— ) BE-condensates, 
and a lower temperature mixed phase (ss) with arbitrary proportions of 2e- and 2h-CPs, see Fig. 9. Of greater 
physical interest are the two higher-Tc pure phases so that we focus below only on them. For each pure phase 
at a critical temperature we have either /S.{Tcs+) ~ /y^io(7cs+) = or lS.{Tcs-) = f \J moiTcs-) = 0, where 
A(T) is the electronic (BCS-like) energy gap. Their intersection gives the BCS Tc value of 7.64 x 10~^Tf that 
also follows from the familiar BCS expression 

Tc/Tf ~ l.lM{hujD/EF)eyi^{-l/\). (68) 

We next focus on s = 1 which, as we saw in Sections 3 to 5, occurs in the leading term for "ordinary" CPs 
in a Fermi sea as well as for "generalized" CPs in a BCS state. For these, the boson energy 77 in (47) to be used 
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later has leading terms in the many-body Bethe-Salpeter (BS) CP dispersion relation that are 

7] ~ {X/27T)hvFK 2D (69) 

[53] and 

r] ~ iX/4)hvFK 3D (70) 

[54]. As before, A = VN{Ep) where N{Ep) is the electron DOS (for one spin) at the Fermi surface. Note that 
the boson energy rj is not the quadratic /2{2m) appropriate for a composite boson of mass 2m moving not 
in the Fermi sea but in vacuum [3]- [6], [11]- [14], and assumed in (48). Note also that the BS CP linear dispersion 
coefficient A/27r in 2D (or A/4 in 3D) contrasts markedly with the coupling-independent 2/7r coefficient in 2D 
(or 1/2 in 3D, as apparently first quoted in Ref. [34], p. 33) obtained [43] in the simple CP problem [29] which 
ignores holes. Thus, again with ub = n/2, for s = 1 and Ci = \b{d)hvF with 6(2) = l/27r and 6(3) = 1/4 
according to (69) and (70), (C.l) leads to 

TJTp = 2Xb{d)/ [dr{d)C{d)f'' . (71) 

For A = 1/2 this is ~ 0.088 if d = 2 since C(2) = and ~ 0.129 if d = 3 since C(3) =i 1.202. These two 

values will appear as the black squares in Figs. 10 to 13 (upper ones in Figs. 10 and 11). They mark the BEC 
limiting values if all electrons in our 2D or 3D many-electron system were imagined paired into noninteracting 
bosons formed with the BCS model interelectron interaction. For A = 1/4 the lower black squares in Figs. 10 
and 11 apply. 

8.1 Two dimensions (2D) 

To address cuprates and other copper-plane-free planar superconductors (such as Sr2Y Rui-xCuxOq with Tc 
49 K, which is Cu-doped but has no CuO planes) we deal first with the CBFM in 2D. Although it is still 
controversial [109] as to which planes, the CuO or BaO or SrO, etc., the superconductivity resides in, these 
planes are parallel to each other, say, in the ab directions. Resistivity p anisotropics pd Pab^ with the c direction 
perpendicular to the ab denoting the CuO or BaO or SrO planes, can be as high as 10"'' in Bi2+xSr2-yCuOf,+5 
[108], if not higher, even though "only" about 10^ in YBa2CuzOT-5 (YBCO). From the Drude 1900 resistivity 
model [110] p = ni/ne^r for current carriers of charge e, effective mass m, number density n, and r is some 
average time between collisions. Thus, if pj Pah = mc/mab is oo one has a precisely 2D situation; if it is 1 
we have the perfectly isotropic 3D case. Hence, the large (10^) but finite ratio observed implies (2 -I- €)D or 
"quasi-2D." For cuprates the value d ~ 2.03 has been extracted independently by two groups [111] [112] for 
YBCO as more realistic than d = 2, since that value reflects inter-CuO (or BaO)-layer couplings. Our results 
in that case would be very close to those for d = 2 since, e.g., from Appendix C (C.l), for s = 1 (but, perhaps 
very significantly, not for s = 2) varies little [113] with d around d = 2. This justifies models in exactly 2D, at 
least as a very good initial approximation. 

In 2D the electronic DOS per unit area is constant, namely N{e) = m/2nh'^. Using the leading term of 
the BS CP linear dispersion relation (69) we get for the bosonic DOS 

M{r]) = {l/2Tr)K{dK/dr)) ~ {2n/X'^n\l)r) (72) 

instead of the constant that follows in 2D from the quadratic dispersion ri = J^K"^ /2{2m). Employing Ef = 
-Kh^nf/m = ksTf as energy/density /temperature scaling factors, and the relation n/uf = {Ep/Ef)'^/'^ to 
convert quantities such as Tc/Tf to Tc/Tp, where Ep = ksTp, the two "working equations" for the pure 2e-CP 
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phase [i.e., niBiTc) = 0] (with all quantities dimensionless, energies in units of Ej and electron particle-densities 
in units of n/) become 



1 + /iwd/S — u = A(/kj£)/2) / dx-. — - — rtanh 



\x - 

2Tr ' 



(73) 



oo 

^2 



dx[l- tanh^] + 3^ - y rfxa;[coth ^^j^^^ ^ - 1] = n. (74) 


These are just the gap-like equation associated with 2e-CPs and its corresponding number equation. Recall 
that Hf can be seen [13] to be the number-density of unpaired electrons at zero T and zero A. 

For the pure 2h-CP phase (i.e., ns = 0) the two working equations are 

1 

/i- 1 -hfej£)/2 = A(/ia;£)/2) f dx- — - — -tanh^-^, (75) 

J \X — ilc 

1 — huJD 

00 00 

- J dx[l - tanh^] -^-J rfxx[coth -^ff^ - ^] = n. (76) 


which are the gap-like equation associated with 2h-CPs and its corresponding number equation. Note that (74) 
and (76) are quadratic in n. Further, all integrals there are exact, i.e.. 



00 

/ 



dx[l - tanh^— ^] =11 + 2Tcln[2cosh(/x/2TJ], (77) 



and 



J dxx[coth^ - 1] = 2T2polyLog[2, e'^/^^] (78) 


where the polylogarithm function as defined in Ref. [114] is 

00 

PolyLog[a,z]=J2^^/l"- (79) 
I =1 

The integrals in (73) and (75), however, were performed numerically. In 2D we use the two extreme values of 
A = 1/4 (lower set of curves in Figs. 10 and 11) and =1/2 (upper set of curves), and hojo/Ep = 0.05 (a typical 
value for cuprates). If A > 1/2 the ionic lattice in 3D becomes unstable, [115] and Ref. [3] p. 204. Equations 
(73) to (76) lead to the Tc/Tp vs. n/uf phase diagram which is graphed in Figs. 10 and 11 for both 2e-CP 
(dashed curve) and 2h-CP (full curve) pure, stable BEC-like phases. The value n/uf = 1 corresponds to perfect 
2e/2h-CP symmetry for the lower-Tc mixed phase (not shown), while n/uf > 1 (and < 1) signifies more (less) 
2e-CPs than 2h-CPs in the mixed phase. The Tc value where both phase-boundary curves no(Tc) = mo(Tc) 
= intersect is marked by the large dots in the figure. The values are consistent with those gotten from (68) 
which gives ~ 0.001 for A = 1/4 and ~ 0.008 for A = 1/2, for hwo/Ep = 0.05. [The values from (68) hardly 
differ from those from the exact BCS (implicit) Tc-formula (Ref. [44], p. 447) 1 = x J^"^^'"''^'' dxx''^ tanhx.] 
Cuprate data empirically [116] show T^s and Tp's falling within the range Tc/Tp ~ 0.03 — 0.09. Thus, moderate 
departures from perfect 2e/2h-CP symmetry enable the CBFM to predict quasi-2D cuprate empirical Tc values, 
without abandoning electron- phonon dynamics — contrary to popular belief. 
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As the variable rif used to scale the electron number density n in Fig. 10 does not appear at all in the 

number equation (57), wc also scaled n with the variable nf{T) evaluated at T = which docs appear in (57) 
and which has a clear-cut physical meaning. Since A(T(.) = 0, from (51) and (52), (60) reduces in 2D to the 
analytic expression 

n/(Te)/n = {T^/Tf) ln[l + exp{/i(Te)/fcsre}]. (80) 

For the pure 2e-CP condensate case, i.e., niBiT) = 0, we must have n/nf{Tc) > 1 from (57), where in general 
71/ (Tc) ^ rif. Here nf{Tc) is the actual number density of unpaired electrons at T = Tc since the number 
equation (57) for mB{T) = then becomes just 

n/(Tc) + 2nB+(Tc) =n. (81) 

Here ns+(Tc) is the total number density of all excited (or K > Q) 2e-CPs, also commonly known as "preformed 
CPs." It is defined in (58) at any T as 

OO 

ns^{T)^jd^M{^)^^^^^. (82) 

0+ 

Preformed CPs seem to be the main conjecture to explain a "pseudogap" [117] observed above in many 
underdoped cuprates. 

The complete number equation (57) then explicitly reads 

2no(T) + 2nB+{T) - 2mo(T) - 2ms+(T) + n/(T) = n (83) 

with mB+{T) being the number of preformed 2h-CPs, defined as the integral term in (59). Using UfiTc) instead 
of n/ to scale electron densities n in (73) to (76) which are solved numerically, scaling all energies with Ep 
instead of Ef as before, and on eliminating n/, gives the results in Fig. 11. Curiously, these curves are hard 
to distinguish from those of Fig. 10. In Fig. 11 the BCS values (black dots in figure) are not associated as 
before with the intersection of the 2e- and 2h-CP condensation curves, but are merely values calculated with 
the BCS Tc formula (68). On the other hand, the Tc values from the CBFM equations scaled with the previous 
(and somewhat opaque) variable n/, match only approximately in the limit n/n/ — > oo with the analytical BEC 
limits of (71). It has been verified numerically that coincidence in the BEC limit is perfect only in the trivial 
limit of zero coupling, as expected. 

In view of Figs. 10 and 11 we note that room temperature superconductivity (RTSC in figures) is possible, 
e.g., in Fig. 11, for n/nf{Tc) sufficiently less than unity. From (83) with nsiT) = no(T) + nB+(T) = this 
translates into this recipe: 2mB+{Tc) must be sufficiently larger than n/(Tc). Specifically, from Fig. 11 RTSC 
occurs for A = ^ if n/nf{Tc) < 0.1, which from (57) means 2mB+{Tc) > 0.9n/(Tc) or that the total number of 
individual holes bound into preformed 2h-CPs at Tc must somehow be larger than about 90% of the number of 
unpaired electrons at Tc. 

8.2 Three dimensions (3D) 

In 3D the electronic DOS if e = H^k^ /2m is the familiar expression 

N{e) = {l/2n^)k''{dk/de) = {mi /2^Tr^h^)^ (84) 

which coincides with (48). Analogously as before, Ef = (?l^/2m)(37r^n/)^/^ = fcsT/ again differs from Ep = 
(?i^/2TO)(37r^n)^/^ = ksTp unless perfect 2c/2h-CP symmetry holds in which case they coincide. On the other 
hand, the leading term in the BS CP boson dispersion energy is now the linear expression (70) so that 

M(r/) = {l/2Tr^)K^{dK/dr]) ~ {32/Tr^X^ti^vl)rj''. (85) 
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The above working equations for the pure 2e-CP phase in 3D (all quantities dimensionless as before in terms 
of Ef and n/) are 



l + n = \{nwD/2)-^ [ dx^i^— tanh^J, (86) 

1 



oo 

^ j dx\/x[l 





while for the pure 2h-CP phase they are 



tanh^] + dxx' [coth (87) 



i 

If 1 

a — 1 + haJD/"^ = MhcJo/^) , ,„ / dx\fx-, rtanh 

n^/'* J \x — ij\ 



\x - ix\ 



. . . ^Tr, 

1 — hoJD 



(88) 



I J dxV^H - tanh^] - / dx.^[coth^i^^ii^£/^ - 1] = n. (89) 



Results in 3D are reported only for A = 1/2, and Hwd/Ep = 0.005 was used here. In Figs. 12 and 13 the 
full curves are again the 2h-CP BEC phase boundaries while the dashed curves are the 2e-CP BEC ones. The 
large dot again marks the BCS T^/Tf values obtainable from (68) of 0.0001 for A = 1/4 and of 0.0008 for 
A = 1/2. Besides Ref. [116], mpirical data for both exotic and conventional, elemental superconductors in 3D 
are also graphed in Ref. [118]. We see that whereas BCS theory roughly reproduces T^/Tp values well for the 
latter SCs, it takes moderate departures from perfect 2e/2h-CP symmetry to access 3D exotic SC T^jTp values, 
which empirically fall within the range 0.01 — 0.1. This is much larger than the range < 0.001 for conventional 
(elemental) superconductors, also shaded in the figure. 



9 Hole superconductivity 

Finally, wc address the unique but mysterious role played by hole charge carriers [46] in the normal state 
of superconductors in general. For example: a) of the cuprates, those that are hole-doped have transition 
temperatures Tc about six times higher than electron-doped ones; and even in conventional superconductors 
[47] b) over 80% of all superconducting elements have positive Hall coefficients (meaning hole charge carriers); 
while c) over 90% of non-superconducting metallic, non-magnetic elements have electron charge carriers. This 
greater "efficiency" of individual, unpaired hole carriers in producing higher TcS is clearly reflected in Figs. 10 
and 11 for both 2D and 3D superconductors, at least insofar as pure 2h-CP BE condensates exhibiting higher 
TcS than those associated with pure 2e-CP BE condensates. 



10 Conclusions 

This review began with a survey of "ordinary" and "generalized" Cooper pairing, and stressed that if hole 
pairing is treated on an equal footing with electron pairing the original "ordinary" Cooper problem (based on 
the pure-kinetic-energy unperturbed Hamiltonian) is meaningless. However, "generalized" Cooper pairs defined 
in terms of the Bethe-Salpeter equation including both electron and hole pairs, in the ladder approximation, 
but based on the BCS ground state as unperturbed Hamiltonian instead of on the ideal Fermi gas sea, gives rise 
to physically meaningful positive-energy resonances with a finite lifetime for CMM ii' > and infinite lifetime 
for if = 0. 
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It was then sketched how five statistical continuum theories of superconductivity — including both the BCS 

and BEC theories — arc contained as special limiting cases within a single theory, the "complete boson-fermion 
model" (CBFM). This model includes, for the first time, both two-electron and two-hole pairs in freely variable 
proportions, along with unpaired electrons. Thus, the BCS and BEC theories are fully unified within the 
CBFM. The BCS condensate [specified by a single equation, namely the T-dependent gap equation for A(T)] 
follows directly from the CBFM as a BE condensate through the condition for phase equilibria when both total 
2c- and 2h-pair number, as well as their condensate, densities arc equal at a given temperature and coupling - 
provided tlie coupling is weak enough so that the electron chemical potential //, rougiily equals the Fermi energy 
Ep- Ordinary BEC theory, on the other hand, is recovered from the CBFM when hole pairs are completely 
neglected, the BF coupling / is made to vanish, and the limit of zero unpaired electrons is taken. 

Lastly, the practical outcome of the BCS-BEC unification via the CBFM is fourfold: a) enhancements in Tc, 
by more than an order-of-magnitude in 2D, and more than two orders-of-magnitude in 3D, are obtained with 
pure electron-pair, and even more so with pure hole-pair, BE condensates for the same clcctron-phonon dynamics 
mimicked by the BCS model interaction; b) these enhancements in fall within empirical ranges for 2D and 
3D "exotic" SCs, whereas BCS Tc values remain low and within the empirical ranges for conventional, elemental 
SCs using standard interaction-parameter values; c) hole-doped SCs are predicted in both 2D and 3D to have 
higher T^s than electron-doped ones, in agreement with observation; and finally that d) room temperature 
superconductivity is possible, with the same interaction parameters used in BCS theory for conventional SCs, 
but only via hole-pair BE condensates. 
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Figure Captions 

1. Dimensionlcss CP excitation energy Sk/{—£o) = ^k/^o i>s K/kp , numerically determined from (8) 
for different couplings B^/Ep, full curves, for the delta interaction in 2D. The dot-dashed line is the linear 
approximation (virtually coincident with the exact curve for B2/EF % 0.1) while the dashed curve is the 
pure quadratic term of (9). Dots denote values of CMM wavenumber where the CP breaks up, i.e., where 

fx =0. 

2. Cross-section of overlap volume in fc-spacc (shading) where the tip of the relative CP wavevector k 
= ^(ki — k2) must point for the attractive BCS model interaction (15) to be nonzero for a CP of CMM 

/lK = fi(ki+k2). 

3. Wavcfunction Feynman diagrams for 2p (■0+), 2h {tp_) and ph (-00) bound states arising from the BS 
equations. Shaded rectangles designate diagrams that do not contribute in the IPG-based case. 

4. a) Exact "moving" 2p-CP (real) energy £k (in units of Ep) in 2D from (37) (full curves), compared with 
its linear leading term (thin short-dashed lines) and its linear plus quadratic expansion (long-dashed curves) 
both from (38), vs. CMM wavenumber K (in units of kp), for interaction (20) parameters A = j (lower set 

of curves) and i (upper set of curves), and hjJo/Ep = 0.05. For reference, leading linear term (34) of trivial 
ABH sound mode is also plotted (lower thick dashed line), b) 2p-CP lifetime as defined in (39). c) Analogy of 
ordinary and BCS-based-BS 2p-CPs with various confined states in a 3D potential problem. 

5. BE condensate- fraction curves 1 - {T/T^y/' for bosons in d = 3, 2, or 1 with dispersion relation (1) with 
s = 2 or 1, for a pure phase of either 2e- or 2h-CPs as discussed in text, compared to empirical data for 3D 
SCs {Nb/Cu and Sn), two quasi-2D SCs (F123 and Bi2212) and a quasi-lD SC (4 A- wide nanotubes). Data 
for the latter are for A(T)/A(0) but are plotted as [A(T)/A(0)]^ so as to reflect 2h-CP condensate fraction 
mo(T)/mo(0) according to (67). The curve marked 1/1 strictly corresponds to Tc = 0; however, it serves as a 
lower bound to all curves with d/s = (1 -|- e)/l for small but nonzero e for which > 0. The ordinate axis is 
labeled with the superelectron number density ns{T) in units of the normal electron density n. Also shown for 
reference are the two-fluid model [80] curve 1 — {T/Tc)^ and the BCS gap A(T)/A(0) order parameter curve 
[81]. 

6. BF vertex diagrams depicting disintegration of 2e-CPs into two unpaired electrons; 2e-CP formation 
process from the latter; formation of a 2h-CP from two unpaired holes; and the disintegration of a 2h-CP into 
the latter. All four processes are contained in Hint Eq. (42). 

7. Flow chart of how the CBFM reduces in special cases to the statistical, continuum models of supercon- 
ductivity discussed in text, thereby displaying how both BCS and BEC theories are unified. The three "legs" 
refer to: perfect 2e-/2h-CP symmetry (middle leg, giving rise to the BCS-Bose crossover and BCS theories); to 
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no 2h-CPs present (right leg) leading to ordinary BEC as a special case; and to no 2e-CPs present (left leg) 
predicting room temperature superconductivity (RTSC), see Figs. 10 and 13 below, via a 2h-CP BE condensate. 

8. Critical 2D BEC-likc temperature (Tc in units of Tp) approached from above within the ideal BF model 
(IBFM) with the BCS model interaction for A = 1/2 for varying fkVD/Ep = Qd/Tf- Results are for: the 
pure unbreakable-boson gas with some and with all fermions paired; for the breakable-boson gas; and for the 
boson- fermion mixture in thermal/ chemical equilibrium (thick full curve labeled "binary gas"), all as described 
in Ref. [18] for original (simple) CPs where Ci = {2/iT)hvF in (1). Dashed curve is the BCS theory Tc, given 
by (68). Cuprate experimental data arc taken from Ref. [116]. The IBFM can be considered a model for the 
normal state, and its temperature instability a prediction of Tc- 

9. Phase diagram [14] with 3D superconducting critical temperature phase boundaries Tcs+i Tcs-, Tcss+i 
and Tcss- as functions of An = n/n/ — 1 as defined in text, in the vicinity of the approximate BCS Tc value, 
assuming the quadratic boson dispersion ry = h?K'^/2{2m), for A = 1/5 and Kwd/Ep = 0.001. 

10. Phase diagram in 2D temperature T (in units of Tp) and electron density n (in units of n/ as defined 

in text) showing the phase boundaries of T^s for pure 2c-CP BEC phases (dashed curves) determined by 
A(Tc) = f^/nf){Tc) = and pure 2h-CP BEC phases (full curves) given by /\{Tc) = f^/mATc) ee for A = 1/4 
and 1/2 with tuvo/Ep = 0.05. Intersections corresponding to no{T) = mo{T) approximately reproduce the 
BCS Tc as given by (67) and are marked by black dots. Black squares mark the BEC limit where all electrons 
are imagined paired into 2e-CP bosons, as calculated in (71), to which values the dashed curves reach only 
approximately in the limit n/n/ oo. 

11. Same as Fig. 10 except that n/nj is eliminated numerically and abscissa refers to n/nf{Tc), where 
nf{Tc) is the actual number of impaired electrons at T,, as given by (80). Here, the intersection corresponding 
to no(r) = mo(r) occurs at T = and not at the approximate BCS Tc values as in Fig. 10. Open circles mark 
values of n/nf{Tc) corresponding to n/uf = oo. 

12. Same as Fig. 10 but in 3D, and only for A = |. 

13. Same as Fig. 11 but in 3D, and only for A = |. 

Appendix A. Why CPs are bosons while BCS pairs are not. 

This assertion is simply visualized qualitatively with the aid of Fig. 2 for the BCS model interaction (15). 
The vector k ends in all points of a simple-cubic lattice in fc-space with lattice spacing 2tt/L with L the system 
size. In 3D these points are within the overlap volume (shaded in the figure) where the interaction is nonzero. 

But, in the thermodynamic limit there arc infinitely many acceptable k values allowed for the energy £k of 
the CP state, whether for ordinary CPs (7), (16) or for Bethe-Salpeter CPs (37). This implies BE statistics for 
either ordinary or BS CPs as each CP energy £k level has no occupation limit. Not so for BCS pairs as we now 
show. 

More quantitatively, consider fermion creation aj^^^ and annihilation a^^s operators that satisiy the anti- 
commutation relations. 



The BCS pair annihilation and creation operators, respectively, are then 

&kK = akaittkiT and ^iLk = '^kit'^Li- 



CA-l) 



(A.2) 



Here 



k = -(ki - ka) and K = ki + ka 



(A.3) 
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are the relative and total (or center-of-mass) momenta wavevectors, respectively, associated with two fermions 
with wavevectors 

ki=K/2 + k and k2=K/2-k. (A.4) 
We show below that 6kK and b^^j^ satisfy: a) the sometimes called "pseudo-boson" commutation relations 

[&kK, ?>k/K] = (1 - "K/2-ki - «K/2+kT)^kk', (A.5) 

[&lcK,''LK] = [&kK,6k'K]=0, (A.6) 

where 

nK/2±kl = aK/2±ki'^K/2±ki (A.7) 

are fermion number operators; as well as b) a "pseudo-fermion" anti-commutation relation 

{fokK, 6k'K} = 26kK6k'K(l - ^kk')- (A.8) 

Our only restriction is that K = ki + k2 = k'^ + k2. Clearly, BCS pairs are not bosons as they do not satisfy 
(Ref. [34] p. 38) the ordinary boson commutation relations 

[^kK, 6iL'k] = '^kk' and [blj^, b^j^] = [6kK, 6k'K] = 0. (A.9) 

If K = (so that ki = — k2 = k), and calling 6kK=o = ^'k, etc., (A.5) and (A.6) become Eqs. (2.11) and 
(2.12) of Ref. [7], and (A.8) becomes Eq. (2.13) thereof. 

To prove (A.5) to (A.8) write 

[bkK, 6k'K] = ^kK^k/K - ''k'K^kK = Ck^ i Clki T Ok'J '^kU ~ 4'J "^i '^ks i Oki T (A.IO) 

using (A. 2). This can alternately be rewritten, using (A.4), as 

[^kK, ^k'K] = aK/2-kiaK/2+kTaK/2+k'T'^K/2-k'i ~ aK/2+k'T'^K/2-k'i"K/2-kiaK/2+kT- (A.ll) 

The pair of Fermi operators (iK/2+kT<^K/2+k'T term gives, using (A.l), 

«K/2-ki(<5kk' - "K/2+k'T"K/2+kT)«L/2-k'i 
= «K/2-kJ.aK/2_k'i^kk' - aK/2-kJ.aK/2+k'T'^K/2+kT«K/2-k'J.- (A-12) 

Using (A.l) again the first term becomes (5kk'(l ^ '^K/2-k'i'^K/2-kj.)- The last term of (A. 12), after anti- 
commuting all creation operators to the left and recalling that S^i = 0, etc., gives — <5kk'C'K/2+k'T'^K/2+kT 
+ '^K/2+k'T'^K/2-k'j.'^K/2-kiaK/2+kT- Inserting this in (A.ll) leaves precisely (A.5) if we recall the number 
operators (A.7). To prove (A.6) write 

[&kK,^k'K] = ^'kK^'k'K — ^'k'K^'kK = ^ka i 'J ki T C^k^ i C^ki T ~ Clk^ i Clk^ T '^ka i Clki T (A. 13) 

using (A. 2); or equivalently 

[bkK.,bk'K] = aK/2-kJ.aK/2+kT^lK/2-k'iaK/2+k'T ~ aK/2-k'iaK/2+k'T'JK/2-kj.aK/2+kT. (A.14) 

The first term is easily brought into a form cancelling the last term by simply anticommuting operators with 
primed subindices to the left, thus proving (A.6). Finally to prove (A.8) write 

{6kK, &k'K} = ^kK^k'K + 6k'K&kK = aksiOkiTak^iOk'iT + ak^i^k'iT'^ksiakiT (A.15) 
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= «K/2-ki«K/2+kT^K/2-k'iaK/2+k't + 'iK/2-k'i'^K/2+k't'^K/2-kJ.<^K/2+kT> (A.16) 

= 2aK/2-kj.aK/2+kT<^K/2-k'i"K/2+k'T = 26kK&k'K if k 7^ k'. (A.17) 

However, if k = k' then 

{&kK,&kK} = 2aK/2-kj.aK/2+kT«K/2-kiaK/2+kT = (A.18) 
since (aK/2±kj.)^ = 0. Hence (A. 8) is true. 

The main point here is simply this: any number of CPs with a definite K but different values of k can 
occupy a state of CP energy £k and thus not only obey BE statistics but (as this in turn demands) also obey 
the boson commutation relations (A. 5) for k ^ k'. Hence, CPs can suffer a BEC as a given CP state involves 
no two BCS pairs with the same k. 

Appendix B. CBFM and BCS condensation energies compared. 

Using (47) for T = when no(T) = mo(T) and ns(T) = ms(T) we obtain for the superconducting state, 
calling A(T = 0) = A, 

O,(T = 0) 



J — u 



where ^ = e — /x. Recalling the expression (52) for A, using (43) and (44), and since /j, = Ef^ gives 

= 2fe^^no(0) + 2 / dmOi + N{Ep) ^^[^-^f + A^]. (B.l) 

J— /J, J —HWD 

The first and second members of the last term have respectively odd and even integrands, so that this term 
becomes 



rhoJD I 

Jo 



'0 

On the other hand, for the normal state no = and mo = 0, so that from (52) A(T) = 0. Hence, from (64) the 
CBFM condensation energy Eg — En per unit volume is just 



Eg — En 



2hiODno{0) + 2N{EF) d^{^-Jf + A^) (CBFM). (B.2) 

Jo 



^From Eq. (2), Sec. 4.3.3.1, Ref. [119] the integral evaluates to 



--hwoVi^^ny + A2 + + -A2 In 



= 1a2 \n{A/2hjD) - - ^[A^/ihjof] + O [AP/{hjnf] ■ 



(CBFM). (B.3) 



Thus, since uq = A^//^, (B.2) for weak coupUng A = ^ becomes 

. - -N(Ef)A^ 

L'^ \^o 2 ^ ' 

By contrast, the original expressions (2.41) and (2.43) in Ref. [7] give 

.hu^o I \ A^ 

= mEril ^[(-^^^y- (BCS) (B.4) 



Es — En 



where V is defined in (15). Multiplying (32) by A^/2 is equivalent to 

rh(Ajo ^2 ^2 



2N{Ef) / d^ 
Jo 



2^/eT^ ^ 
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which when combined with (B.4) gives 



Eg — En 



= 2N{E^) / 



huJr 



1 2r + A2 



Using Eqs. (3) and (10), Sec. 4.3.3.1, of Ref. [119] finally gives 



Eg — E„ 



= N{EF){nwDf 



which on expansion leads to 

Ec, — En. 



N{Ef)A'^ 



1 - - {A/huJof + 0{A/fwjDf 



(BCS). 



(B.5) 



(B.6) 



Thus, the CBFM condensation energy (B.3) is larger than that of BCS, but only the latter is a rigorous 
upper bound to the exact condensation energy since it follows [7] from a variational trial wavefunction. 

Appendix C. BEC as limit of all electrons paired. 

Here we discuss generalizations of the well-known result that the BEC transition temperature of a d- 
dimensional A/'-fermion gas of Fermi temperature Tp in which all fermions are imagined paired into bosons, is 
just 0.218Tp. (See dashed line in "Uemura plot" of Ref. [118], Fig. 2). These results will provide a convenient 
pure BEC limit which in effect turn out to be an upper limit for the pure 2e-CP phase separation boundary Tc 
values deduced in Sec. 8. 

The general BEC Tc-formula for identical noninteracting bosons in d-dimensions of energy rj = Cs K^, s > 0, 
where as before K is the boson CMM, is [16] 



9± 



sT{d/2){2-KY 



s/d 



where ub is the boson number-density and ga{z) the Bosc integral [120], z = e'^'^/^'^'^ is the "fugacity^^ and 
fj,Q the boson chemical potential. For z — 1, 5<t(1) = C(''')i the Ricmann Zeta-function, if cr > 1, while for 
< a < 1 the infinite series 5cr(l) diverges. Eq. (C.l) is formally valid for all d > and s > 0. Hence, for 
< d < s, Tc = since ffd/s(l) = oo for d/s < 1 but Tc is otherwise finite. We stress that as a consequence of 
the former all 2D Tc predictions in Fig. 7 (except the BCS one that survives for all d > 0, including d = 1) 
would collapse to zero had s = 2 been used in 2D instead of the correct s = 1 arising from the Fermi sea. For 
s = 2, C2 = h'^/2mB (C.l) leads to the familiar 3D result 



Tc ~ SMfi^n^J^ /mskB 



(C.2) 



since C(3/2) ~ 2.612. Recalling that fesTp- = h^k%/2m with kp = [2'^-'^-K'^/'^dT{d/2)nY/'^ from (36), then if 
niB = 2m and ub = n/2 (all electrons paired) for s = 2 (C.l) implies that 



Tc/Tf = i[2/dr(d/2)5rf/2(l)]2/'' =0 for d<2 



(C.3) 



since 5^/2(1) = 00 for d/2 < 1. When d > 2 then Tc is nonzero. For d = 3 we get the familiar limit mentioned 
before 

Tc/Tf = \ [2/3r(3/2)C(3/2)]^/^ ~ 0.218, (C.4) 



(see also dashed line Fig. 2 of Ref. [118]). 
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